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Abstract 

We construct time dependent random fields on the sphere through coordinates 
change and subordination and we study the associated angular power spectrum. Some 
of this random fields arise naturally as solutions of partial differential equations with 
random initial condition represented by a Gaussian random field. 

Keywords : random field, power spectrum, coordinates change, time change, subordi- 
nator, fractional operator. 

1 Introduction and main results 

In this paper we construct random fields on the unit sphere with different covariance struc- 
tures by means of the coordinates change of Gaussian random fields given by the subordinate 
Brownian motions on the sphere. We find out random fields with spectral representations 
on the unit sphere which are solutions to partial differential equations, with random initial 
condition, involving the fractional operator 



ni fix) = / (/(y) - fix)) J(x, y)\idy) 

where A is the Lebesgue measure on the unit sphere . Such random fields turn out to be 
interesting because of their covariance structure and, in particular they well explain certain 
forms of the angular power spectrum studied in literature. In our view, the present work 
directs one's efforts towards the analysis of the Cosmic Microwave Background (CMB) radi- 
ation which is currently at the core of physical and cosmological research ( [31 [9l [15] ) . Indeed, 
our model well accords with the results in [21 \TT\ \T2\ [T3] where the authors deal with high 
frequency (or high resolution) analysis of random fields and, in particular with [13] where 
a special class of models depending on the angular power spectrum has been introduced 
and relationship between ergodicity and asymptotic Gaussianity has been investigated. The 
mathematical objects we introduce here are time dependent random fields which well depict 
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the random movement of particles in random environments. In particular, if T{'d, (p) is a 
random field on the two dimensional sphere, then as time passes, 

{^t,iPt,T{^t,ipt)) en^ t>o 

is a randomly varying point in a random environment given by the stochastic sphere {i^t, ft), 
t > 0, which is the rotational Brownian motion or Brownian motion on the 2-sphere. Since 
the processes we deal with are random solutions of Cauchy problems, we also know their 
evolution in time. The angular power spectrum of the random fields considered in this work 
exhibits polynomial and/or exponential behaviour in the high frequency (or resolution) 
analysis and therefore, we introduce a large class of models in which many aspects can be 
captured, such as Sachs- Wolfe effect or Silk damping effect for instance. 

1.1 Statement of results 
1.1.1 Preliminaries 

We focus on the zero-mean, isotropic Gaussian random fields T on the sphere with spectral 
representation 

OD +1 OO 

/=0 m=-l 1=0 

where 

aim= [ T{x)Yi*^{x)X{dx) (1.2) 

JS2 

are Fourier random coefficients. We consider the square integrable n- weakly isotropic (with 
rotationally invariant law) random field {T{x); x G Sf } on the sphere = {x E : \x\ = 
1} for which 

ET{gx)=0, and T{gx)^=T{x), for ah 5 G 50(3) 

where represents equality in law of stochastic processes and SO{3) the special group 
of rotations in R^. The random field T is said to be n-weakly isotropic if K\T(gx)\^ < 00 
(n > 2) for every x G Sf and if, for every xi, . . . , x„ G Sf and every g G SO{3) we have that 

E[T(xi) X • • • X T{xn)] = nrigxi) x • • • x T{gxn)]. 

Convergence in (jl.ip holds in the mean square sense, both with respect to L?'{dP (8> \{dx)) 
and with respect to L'^{dP) for fixed x £ Sf, X{dx) is the Lebesgue measure on the unit 
sphere Sf and the set of spherical harmonics {Yim '■ I > 0, m = —I, . . . , +1} represents an 
orthogonal basis for the space -L^(Sf , A((ij;)), with X{dx) = X{d'd,dip) = dipd'd sini? being 
X G Sf represented as 

X = (sin(^cos'!9,sini?cos(/?,cosi?), G [0, vr], 93 G [0,27r]. 

Sometimes, we will also write f{x) = f{'d,(p). The random coefficients (|1.2|) are zero-mean 
Gaussian complex random variables such that ([1]) 

E[ai^a^,^,] = SUtCi = Ejaz^p (1.3) 
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where C/, / > is the angular power spectrum of the random field T which fully 
characterizes the dependence structure of T (under Gaussianity) and 

^6 ^ f 1' a = & n 4^ 

" \ 0, a^b ^ ' 

is the Kronecker's delta symbol. The symbol "*" stands for complex conjugation. The 
interested reader can find a deep discussion and presentation of results concerning this field 
in the book by Marinucci and Peccati |14j . 

For a fixed integer I and m = /(/ + !), the spherical harmonics 



(or linear combination of them) solve the eigenvalue problem 

As2Yim = -tJ-iYim, l>0,\m\<l (1.5) 

where 

is the spherical Laplace operator. We also recall that 



are the associated Legendre functions and the Rodrigues' formula 

1 



defines the Legendre polynomials. 

We introduce the transition function 

oo +1 

^tfix - xo) = E^fiBf -xo) = J2Yl lzm(x)y4(^o)Pf (w) (1.7) 

1=0 m=-l 

where 

P*(/i) = exp i-t^i^i)) = Eexp (-^A) (1-8) 

is the Laplace transform of the law of the subordinator Dt, t > with Laplace exponent ^. 
Formula ()1.7p with ^{^) = fJ, (and therefore Dt = t, the elementary subordinator) is the 
solution of the Cauchy problem 

(S.-A,;)„ = 

uo = j{x- Xo) 
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and, for f = 6 (the Dirac delta function, this means fim = 1 for ah l,m) is the transition 
law 

°° 2/ + 1 

u{x, xo,t) = Euo{x + Bt-xo) = Y^ -J^Qiii^^ xo))e-*'^' (1.10) 

1=0 







of the rotational Brownian motion Bt, t > 0, with values in Sf and starting point B{ 
xo € Sf . In formula (jl.lOp we have used the addition formula (|2.16p below where 

{x,y) = cos d{x,y) 

is the inner product in R'^ and d{-, •) is the usual spherical distance. Thus, formula (|1.7p is 
the spherical convolution of / with the transition law (jl.lOp of the subordinate rotational 
Brownian motion Bf = Bjj^, t > 0, with Bq = xq. 

Prom now on, we will consider the subordinate process Bf — xq with law ()1.7p and 
flm = 1 for all l,m. 

Assumption 1. For all ti,t2 > and x,y G Sf, the subordinate rotational Brownian 
motions Bf^ — x and Bf^ — y are independent if and only if x ^ y. 

Remark 1. We notice that this assumption is rational and justified by the fact that for dif- 
ferent starting points, the processes Bf^ and Bf^ are two independent copies of Bf whereas, 
if X = y, then this means that i?* and i?* are two realizations of Bf at different times ti 
and t2 ■ 

1.1.2 Subordinate stochastic sphere 

The main object we deal with in this work is the time-dependent random field 

Tf(x) =T(B* -a;), t > 0, x G S? (1.11) 

which is the spherical random field T indexed by the subordinate rotational Brownian 
motion Bf , t > 0, starting form x G Sf. We say that the set {Bf , t > 0} C Sf is a 
subordinate stochastic sphere and represents the new set of indices by means of which the 
random field T is observed. 

Let 5^5* be the cr-field generated by Bf . We point out that ([6]) 

E[Tf(x)]" =-E[¥.[T{Bf - x)\^sA = ^ \^[T {gx)r\^ b^] = WigxT (1-12) 
for all g € SO{2>) due to the fact that ([H]) 




(which does not depend on x) where E[a;^o ' ' ' (^i„o\ is called the angular polyspectrum of 
order n — 1 associated with the field T. We also recall that 

E[r(x)]2 = ^^Q<c^ (1.13) 

1=0 

and T is square integrable. 
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Theorem 1. (Space-Time covariance) For < ti < t2 < oo and x,y £ Sf, we have 
that 

oo , 

E[T,^(x)T*(y)] =J2^QQi{{x,y))P,%^{f,0, x^y (1.14) 

1=0 

where Ci is the angular power spectrum of the random field T. 

We recall that the angular power spectrum of T is given by (jl.3p from which we recover 
the covariance 

oo +1 °° o; _|_ 1 

nT{x)T{y)] =Y,Ci YUx)YL{y)=Y.^CM{x,y)). (1.15) 

1=0 m=-l 1=0 

In the last step we have considered once again the addition formula (|2.16p . Let us write the 
space covariance ()1.14p as 

It = lt{x,y) = E[ll{x)%l{y)], for x ^ y, < t = ti + t2 < oo. (1.16) 

We can immediately verify that 

\im jt(.x,y) = 7o(x,y) = K{T{x)T{y)] 
t^o+ 

and 

C-o 

lim 7t(x,y) = 7oo = hm 7t(x,y) = — . 
t— s-oo t-^-oo 47r 

Indeed, 7oo immediately follows from (jl.Sp whereas, from the representation (jl.ip and for- 
mula (jl.3p we get (jl.lSp . Therefore, the time-dependent random field (jl.lip is the Gaussian 
random field T at time zero and becomes an uniformly distributed random field as time 
passes and goes to infinity. 

Theorem 2. (Time covariance) For < ti < t2 < oo and x,y G Sf, we have that 

E[TZ{gx)%Z{gx)] = ^ (w), V<7 € SO{3). (1.17) 

1=0 

Remark 2. The variance of%f {x) can he obtained from (|1.17p hy observing that Pq {ni) = 1 
for all I > 0. Also we observe that for the covariance 

7, = 7,(x,x) =E[T*(x)Tt*(x)], 0<r = t2-ii<oo (1.18) 

we get that 700 = Co/47r. 

Remark 3. We notice that < (fii) < 1. 
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We say that the zero- mean process Xt, t >0, exhibits long range dependence if 

oo 

T=l 

Otherwise, we have short range dependence. In the following remarks we focus on the high 
frequency behaviour and the dependence structure of the random field considered so far. 
Since 

oo 
1=0 

the random field 

+1 

'^fA^) = E "'imyUBf -x), X G S?, t > 0, / > (1.19) 

represents the projection of on T-Li for each / > and, the high resolution analysis 

(for large I) of T*(x) reveals interesting properties of T*(x). As I increases we get more and 
more information about observations. We say that ()1.19p is the high frequency component 
of T*(a;). Furthermore, (see formula ()2.13p below) 

oo oo 

E[^lix)lZ{x)] = E[ll^ix)lf,^ix)] = J2mf,tA^)'^tA^)] (1-20) 

l,l'=0 1=0 

which means that (jl.lOp are uncorrelated over /. 

Remark 4. Let ^ be the symbol of a stable subordinator, ^'(/u) = We have that 

f;E[T;^,^,(x)T5(x)] =^-^C,f^e~^^r 

r=l r=l 

2/ + 1 ^ 1 



Att e^i - 1 

21+1 _,2a 

for large I. We have the same behaviour if we consider a stable subordinator with drift, that 
is ^(n) = bfj. + fj,°' with b> 0. 

Remark 5. Let ^ be the symbol of a geometric stable subordinator, ^(^) = ln(l + yU°). We 
have that 

oo 



9/ _L 1 1 

r=l r=l ^ ' n ' 



_2/ + l^ 1 

for large I. For a = 1, the geometric stable subordinator becomes the gamma subordinator. 
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Remark 6. The power spectrum is commonly assumed to be Ci ^ l""' for 7 > 2 to ensure 
summability (in (jl.lSp for instance). Our power spectrum, say Ci, takes the form Cie~^ " 
or Cil~'^°' as shown in the previous remarks and therefore, introduces a class of models in 
terms of power spectrum of Gaussian field ( a slightly modification of the class D introduced 
in 113^ ) in which, for G R, € [0, 00), 

00 

J^^Z-^e"'" < 00, L large. (1.21) 

l=L 

We notice that v and 6 must be such that if u = {), then 9 > 2. In U3^ . the authors have 
considered the class D in which they found connections between high frequency ergodicity 
and high frequency Gaussianity (for large I ). 

Remark 7. Under the assumption that Ci ~ /"'^, 7 > 2, the random fields considered in 
Remark [7] and Remark exhibit short range dependence but the rate of convergence of the 
sum over I is different. This means that the I th frequency component of1^f{x) has different 
covariance structure as /—)■ 00 for different symbol ^. 

Remark 8. Let us consider the sum of an a-stable subordinator, say Xt and, a geometric 
j3-stable subordinator, say Yt. We assume that Xt and Yt are independent. For c,d>0, the 
subordinator X^ + Y^t, t > has the symbol 

^(^) = c/i" + (iln(l + /), a,/3G(0,l). 

The corresponding covariance function of%f{x) leads to 

°° 9/ _L 1 °° / p-c^l? \ ^ 



4^ ^V(i+/^f)' 
21+1^ 1 

-U; 



4vr e^/^r + Hi^e'^^'T - 1 
21 + 1 2dl3-cP°' 

^^Q/ e . 

We immediately recognize the class we dealt with in Remarkl^ 

In [6] we have considered the time-changed Brownian motion on the sphere -B^j-, t > 0, 
ly G (0, 1) as coordinates change for the random field T. The inverse process S,^ can be 
regarded as the hitting time of a stable subordinator, say Sj^, and therefore written as 

mm) 

£r = inf{s>0 : i^^ ^(0,t)}. 

The resulting field is related to equations involving fractional derivatives in time of order 
v and, it has been shown that, its covariance structure reveals long range dependence for 
the high frequency components. In particular, the covariance function of the Ith frequency 
component is written as 



for large I, where Ei, is the Mittag-Leffler function. As we can immediately see, the sum 
over T of the covariance above is infinite because of z/ G (0, 1). 
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1.1.3 Mean subordinate stochastic sphere 

We now introduce the zero-mean Gaussian random field 

r?f(x)=E[Tf(x)|5T] (1.22) 

where is the o"- field generated by the random field T on Sf. We say that the random field 
T is indexed by a mean subordinate stochastic sphere meaning that (jl.22p can be written 
as (see formula (|2.15p below) 

oo 

r]fi^)=Yl E aimmmiBf-x)] 

1=0 \m\<l 
oo 

=E E «/™e-*'^(«)i^^(x) 

^=0 |m|<« 

oo 

= ^e-ti'(w)7|(x), xeSlt>0. 

1=0 

Furthermore, the covariance function is written as (see formula ()2.19p below) 

1=0 ^ 

which coincides with p.l4p . Since the covariance above is defined for all x, y, we can write, 
for X = y, 

E[r,figx)f = Y,^Qe-''^^>^^\ V5 G 50(3) (1.23) 
1=0 ^ 

which does not coincide with E[T*(a;)]^ in (jl.lSp as it is in (jl.l2p for T*(x). The variance 
(|1.23p can be also obtained from the higher-order formula 

nvfia-r = E ^""^"^^J^^^^nai.o---auo] (1.24) 

for all g € 50(3) and n € N, which also differs from the higher-order moments of T. As we 
can see, in this case Assumption [1] fails and the random field ()1.22p is qualitatively different 
from (jl.lip . Since ^' > (^' is a Berstein function), we notice that E[ry^(a;)]" ^ as t cx) 
for all n. The Fourier random coefficients are obviously time-dependent and such that 

which comes directly from ()1.2p . We can also write 

E\aim{t)\^ = e-^^'^^'^'^Q, I > 0, m = -I, . . . ,+l. (1.25) 
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From the Parseval's identity we get that 

„ oo +1 oo +1 

•^^1 1=0 m=-l 1=0 m=-l 

or equivalently we say that rjf {x) € L'^{dX (g) dP). 

Let / be a square integrable function on the unit sphere, / € L^(S^). For a € (0, 1), we 
introduce the fractional operator 



/•OD 

fix) = / (P, fix) - fix)) Mids) (1.26) 
Jo 



where M(-) is the Levy measure introduced in Section [2] below and Pt = e is the 
transition semigroup of the rotational Brownian motion Bt- As we will show below, the 
operator (|1.26p can be rewritten as 



nifix)= ifiy)-fix))Jix,y)Xidy) (1.27) 



where A is the Lebesgue measure on and 

with ^'(/i) = d/dfi^ifi). 

The next results are devoted to the study of the solutions of certain partial differential 
equations with random initial condition. Seminal works in this field are those by Krav- 
varitis [H] and Kampe de Feriet [7]. The literature concerning random solutions of partial 
differential equations is vaste and therefore we skip a further introduction and jump to the 
following results. 

Theorem 3. Let ^ be the symbol of a subordinator with no drift. For a € (0, 1), the 
solution to 

idt-m,,)r^fix)=Q 
Vo = T{x) 



(1.28) 



is a Gaussian random field on S\ which can be written as 

oo +1 oo 

vfi^) = E E ^ummix) = Y,^~'''^''^Tiix) (1.29) 

1=0 m=-l 1=0 

where 

aim{t) = aim exp (-il'(w)) , * > (1.30) 
and the series converges in the sense that for all t >0 

2 

''Ixj 



lim E 

L— ^oo 



(1.31) 
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Furthermore, 

vfix - xo) = ¥tT{x - xo) = E[lf{x - xo)|^t] (1.32) 
where 'St is the a-field generated by T. 

Remark 9. Let ^ he the symbol of a subordinator with drift b > 0. The problem (jl.28p 
becomes 

r idt-bAs2 -Ol,)v?ix) = 

I v!ix) = nx) ^^-''^ 

and the random solution is 

oo 

rjf{x) = e-**(^')Ti(x). (1.34) 

1=0 

where ^ is written as in formula ()2.4p below. 

We now study the special case ^'(^) = that is to consider a subordinate rotational 
Brownian motion with a random time which is a stable subordinator. The corresponding 
random field T^(x) solves a stochastic fractional equation involving the fractional spherical 
Laplacian with random initial condition given by the random field T{x), x € Sf . 



Corollary 1. For a € (0,1), the solution to 



idt + i-As2r)r^tix) 
r]o{x) = T{x) 



(1.35) 



is a Gaussian random field on Sf which can be written as 

oo oo 

rjtix) = ^ e-*^r ^ ai^YUx) = ^ e-*^"T;(x) (1.36) 

1=0 \m\<l 1=0 

in the sense of ()1.3ip . Furthermore, 

r]t{x - xo) = FtT{x - xo) = E[Tf (x - xo)|^t] (1-37) 
where St is the a-field generated by T and ^(^) = 

The definition of fractional powers of operators is given below and intended in the sense 
of Bochner's subordination rule (see for example [TO]). This is an interesting topic with a 
long history, we refer to [5] and especially to the references therein for alternative definitions. 

Remark 10. Let the symbol ^ be as in Remark\^ For a fixed time i > 0, say t = 1 for the 

sake of simplicity, the coefficients (|1.25p define an angular power spectrum, say Ci, which 
takes the form 

Ci = Ci e-(-/^"+'^i-(i+M^)) ^ Ci l-^^^e-^"" (1.38) 

for large I. In particular, we obtain a large class of random fields with angular power 
spectrum (|1.38p and parameters f3,a G [0, 1] such that 

oo 

Ct < oo, L large. 

l=L 
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In the analysis of the CMB radiation for instance, in modelling the Sachs- Wolfe effect we 
have a polynomial decay of the angular power spectrum whereas, the Silk damping effect 
entails an exponential decay as I increases. 

Corollary 2. Let us consider ^(^) = /u", a G (0, 1) in (jl.l6p . We have that 

{dt + (-As2)")7i(x, .) = = (a* + (-As2)")7i(-,2/) (1.39) 
in both space variables x and y, subject to the initial condition 

7o(x,y)=E[r(x)r(y)] 
where x,y G Sl with no restriction. 

2 Auxiliary results and proofs 

Let us consider the Levy process Fj, t > 0, with associated Feller semigroup Pt /(x) = 
E/(x + Ft) solving dfU = Au with initial datum uq = f . For the infinitesimal generator A 
of Ft, t > the following representation holds 

M/)(x) = - e-'^^-^mmd^ (2.1) 

for all functions in the domain 

D{A) = |/ G L2(R^dx) : j^m\m?d^ < oo} (2.2) 

where /(^) = J-^d e*^ ''/(x)dx is the Fourier transform of /, $(.) is continuous and negative 
definite. We say that Pt is a pseudo-differential operator with Fourier symbol Pt = exp(— 1$) 
and, $ is the Fourier symbol of A, that is {Af){^) = —^{^)f{^)- Levy processes with 

HO = + (e^«^ - l) M{dy) (2.3) 

where b > and the Levy measure M(.) satisfies J^{y A l)M{dy) < oo and M(-cx),0) = 
0, are the so-called subordinators and possesses non-negative increments, that is non- 
decreasing paths. The Laplace transform of the law of a subordinator leads to the Laplace 
exponent 

/■oo 

^r(^) = -$(i/i) = bfi+ (1 - e-f'y) M{dy) (2.4) 

JO 

for each > and therefore, we can write P^ {n) = exp (— 1^'(^)). We write below some 
explicit forms of the Laplace exponent ^ and the measure M of the corresponding subor- 
dinator (a G (0, 1)): 

• ^(^) = stable, M{y) = ay~°'~''/T(l - a); 
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• ^'(yu) = bfi + stable with drift, M(-) as above and 6 > in ([2^ : 

• ^'(/u) = ln(l + /i): gamma, M{y) = y^^e^"^; 

• ^'(/u) = ln(l + fi°'): geometric stable, M{y) = ay~^Ea{—y) where 

oo 

^ , s \ - 2:-' 

is the Mittag-Leffler fmiction. 

Let i^f , t > be a stable subordinator with a G (0, 1). The generator of F^^q, t > is 
given by the beautiful formula 

for all / G =5^ (the space of rapidly decaying C°° functions) and we formally write the Feller 
semigroup of Fj, t > as Pj = e*"^ which must be understood in the sense of functional 
calculus. For ^ = — is the translation semigroup (or shift operator in the context of 

operational calculus) and formula ()2.5p becomes the Riemann-Liouville fractional derivative 
governing a-stable subordinators. It is well known that, for $(^) = |^|", formula (|2.ip gives 
us the fractional power of the Laplace operator which can be also expressed as 

-(-AR.)"/(x)=C,(a)p.v. / ^^"'Ifl^/^^'^ dy, x G R"^ (2.6) 

where "p. v." stands for the "principal value" of the singular integrals above near the origin 
and Crf(a) is a normalizing constant depending on d and a. The subordinate Brownian 
motion on the sphere is obtained by considering the rotational Brownian motion on Sf with 
a randomly varying time given by i^f . The governing operator is the fractional power of 
the spherical Laplace operator (|1.6p . that is — (— Ag2)", a G (0,1) or Ps = exp(sAg2) in 
formula ()2.5p . 

Hereafter, we use also the symbols 

00 +1 

Im. I m I \m\<l 1=0 m=-l 

in order to streamline the notation as much as possible. 

First we discuss on the solution to (jl.9p . Let / G L^(Sf). Since {Yim : Z G N, m G 
Z and \m\ < 1} is dense in L^(Sf), the function / (see for example the Peter-Weyl repre- 
sentation theorem on the sphere in [l3] or, the references therein) can be written as 

fix) = f{^,^) = Y,flmYU^,y^) 
Im 
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which holds in the sense, where 

hm= f f{^,v)Y^U^,^) sin^d^dif. (2.7) 

The angular power spectrum of / is defined as 

Ai = Mf)= \fimf. (2.8) 

\m\<l 

By fohowing the same arguments we get 

fix -Xo) = J2 flmyim{x)Yil{xo). (2.9) 
Im 

The strongly continuous semigroup Pf associated to Ag2 is written as 

Pt fix - xo) =E^7(5t -^o) = Yl fimYim{x)Yil,{xo) exp {-t^il) 

Im 

where xq G is the starting point at = of a rotational Brownian motion Bt, t > (or 
Brownian motion on the unit sphere Sf ) and comes out from 

E^f{Bt-xo) = J] j;/,„i1w(x)i^:,(xo) I Yil^,{y)YUy)X{dy) exp(-t//0 

I'm' Im -^^i 

by considering ()2.9p and the orthogonality property of spherical harmonics 

/ YUx) Yir^,{x) X{dx) = 6l' (2.10) 

where Y^^ = {—l)"^Yi^m is the complex conjugation of Yim- For the sake of clarity we notice 
that, according to (fTT]) . Pt = Ft for ^{fi) = fi. 

Remark 11. We observe that: 

1. Po = Id, 

2. = Ftlsi = I, 

3. PtP, = Pt+,. 

Indeed, (1) is easy to check, see (|2.9p ." (2) follows from ()2.10p .' (3) is a consequence of the 
Bochner subordination rule and the fact that the law of Bt satisfies the Chapman- Kolmogorov 
equation. Set fim = 1 for all l,m, that is f = 6 and Ft6 is the transition function of Bf , 
t > 0. For all x,y £ (reproducing kernel) 

[ Qi{{x,z))Qi,{{z,y))X{dz)=5l'-^Qi{{x,y)) (2.11) 

and Ft6{x — xq) satisfies the Chapman- Kolmogorov equation. 
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The convergence of the series expansion (jl.ip is understood in the following sense 

2 



lim E 

L— S-oo 



/ T{x) - ^ ^ aimYim{x) X{da 

•^^^ V 1=0 m=-l I 







and also 



^irn^E {t[x) Y.Y. a'^^'-C^)) = 

V Z=0 m=-l J 

as pointed out in the introduction. The coordinates change random field is defined as follows 

(^) = - x), X G S?, t > (2.12) 



Im 



where the subordinate rotational Brownian motion is a measurable map from (17, 5^^*, P) 
to (Sf , =^(Sf ), A). Therefore, the convergence of ()2.12p must be understood in the sense 
that for all t > 



hm E 

L— >oo 



and also, Vt 



E 



/ (nBf -x)-Y,Y^ aimYUBf - x)] 

■^^^ V 1=0 m=--l J 



\{dx] 



lim E 

L—>-oo 



E{T{Bf-x) 



L +1 \ ■ 

aimYim{Bf - x) 

1=0 m=-l J 



0. 



Prom (jl.ip we write 



where 



Tf(x) = 5]T5(. 



i>0 



"^fM = ^ImYlmiBf - X), t > 0, Z > 



m|<i 

is a random field representing the /th frequency component of or the projection of T* 
on which is the space generated by {YJm '■ |"^| ^ 0- We have that 

E[T*,^(x)T;^,,,(y)] = E EKm.«U]IE[y,,^,(S,^-x)y,:„^(S*-y)]. 

|mi|<ii |m2|<;2 

and the Fourier random coefficients aim are independent Gaussian complex random variables 
for which (jl.3p holds. Thus, we get that 

mttM'ittM] =^\A. E cinyuBi - x)YCm{Bi - y)] 

\m\<l 
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and therefore, the covariance function of can be written as fohows 

E[Tf,(x)T*(y)]=j;E[T5^(x)T5,(y)], x^y (2.13) 

and we recall that Ci, Z > is the angular power spectrum of the random field T. 

Proof of Theorem[J\ We have that x ^ y, then for the Assumption [H we write (j2.13p as 
follows 

mitA^)^Z,tM] E QE[Yi^{Br^ - x)]e[Yi:^sbz - y)] 

\m\<l 

where, from the fact that 



2^ + 1 



Qi'{{x,y))YUy)X{dy) = 6l y,:;,(x), (2.14) 



2/' -I- 1 /■ 

E[y,„« - x)] =J2^- / , yim{y)Qi'{{y,x))X{dy)P,^{f,i,) = ycmp?M) (2-15) 

7/ J Si 



we obtain that 

2/' + l 

/' 

which is as to write fim = '^f'^m' i'^ (|2.7p with / as in formula ()1.7p . From the fact that 
= (-l)™il-m, we also obtain 

nyLiBi - y)] ={-^rnyi~n.{Bi - y)] 
=(-i)-itm(x)p|(w) 

=YUy)P?M)- 

Thus, we get that 

\m\<l 

and, from the addition formula 

2/ _(_ 1 

YUx)YrrrXy) = -^Qi{{x,y)) (2.16) 

|mj<Z 

we arrive at 

2/ -I- 1 --^^ 

E[T,tt,(:r)T*,,(y)]=44^QP,'f+,^(w)QK(^,?/)) (2.17) 

which is the (space-time) covariance of the /th frequency component of T*. By combining 
(|2T7l) with (|2l^ we obtain the claim. □ 
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Proof of Theorem [H We start once again from the fact that 

\m\<l 

as we have shown in the previous proof. Here we have that x = y and therefore, from 
Assumption [H we get that 

21' + 1 



An 



Yim{y2)Yim{y2)Kdy2)Pu^u (w) 



Qi'{{yi,y2))Yim{yi)Yim{y2)Kdyi)X{dy2)P^_t^{ni' 



Thus, formula (|2.13|) takes the form 



E[T,^ (x) %l (x)] = E 



and this concludes the proof. □ 
We now consider the random field on the mean subordinate stochastic sphere 

7?f (x) = E[lf (x)|5t] = WiBf - x)\^t], X G S?, t > (2.18) 

where T is the Gaussian random field and 5t is the ci-field generated by T. From 

(jl.ip and ()1.3p we can write the covariance function of ()2.18p as 



nrili^)ril{y)] =E 



Irn I'm' 



.,Y^Cie-^t.+t.m,dYUx)YC^{y) 



Im 



I \m\<l 



(2.19) 



for all x, y € Sf and ti,t2 > 0. 

Proposition 1. For n G N, the higher-order moments of ()2.18p are given by 
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Proof. The higher-order moments of ()2.18p can be obtained as follows 



where 



E 



The random field T is isotropic. Due to the fact that Ti{x) = T{xi\f) where xn = (0,0) 
is the North Pole and that Yim.{xN) = for m 7^ and Yiq{xn) = v^(2/ + l)/47r f[T7]). we 
can write 



E 



'n-=i(2/. + i) 



(4 



By collecting all pieces together we get 



.^(,,)^/ n"=i(2/. + i) 

V (47r)" 



(2.20) 



By observing that 



E[r^f{9xr=E[rjf{x)r, G 50(3) 
we complete the proof. 

The fractional operator (|1.26|) can be rewritten as 

roo 

nifix)= / {Psf{x)-f{x))M{ds) 
Jo 

{Ef{x + Bs)-f{x))M{ds) 
E [{fix + Bs) - fix))] Mids) 

ifiy) - fix)) Xiy)Pr{x + B,e dyjMids) 



□ 



Jsj 

ifiy)- fix)) Jix,y)Xidy) 
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where 

POO 

J{x,y) = / Pr{x + Bs e dy}/dyM{ds) 
Jo 

o/ _|_ 1 rco 
2/ _|_ 1 

and ^'(fJ-) = d/d^^{^) as we can see from (|2.4|) with & = 0. We notice that 

< J{x,y) < oo and J{x,y) = J{y,x). 

We also recall that ^' is a Bernstein function and therefore ^ > (in particular, ^(0) = 0) 
and {-lfd^/dixH{n) < for all > and A; = 1, 2, . . .. 

Proof of Theorem O First we show that 

¥sYUx) = YU^)e-''''^^^\ (2.21) 

Indeed, we have that 

^tYUx) =^YUx + Bf) 

= Y,PM')Yv^'{^) [ YUy)Yv^.{y)\{dy) 



I'm' 



=p^{f^i){-irYr^^{x) 

=P^{^ii)YUx) 

were we have used the orthogonality property (j2.10p and the fact that Y/m = {—^)^Yj*_^. 
We recall that _ 

P*(;u) = Eexp (-^S^J = exp (-i^(/i)) 

where Dt is a subordinator with Laplace exponent ^ . For ^(/i) = /i, we get that 

PsYi^{x) = Yim{x)e-'''^ (2.22) 

2 

where Pt = e is the transition semigroup of the rotational Brownian motion Bt, t > 0, 
with transition law (jl.lOp . We assume that ()1.29p holds true. We immediately see that 
r]Q{x) = T{x) which is in accord with the initial condition. With (|2.22p in mind, we obtain 



mvfix) = / {Ps rjfix) - r?f (rr)) M{ds) 
Jo 

/•oo 

= V e-*"^(''')a;™ / {Ps Yi^{x) - Yi^ix)) M{ds) 
Tr ^0 
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poo 

= y^e-'^^^^^airn. / [e-'^'Yim{x) - YUx)) M{ds) 

POO 

= Ve-**('^')ai^ll™(x) / (e-^^' - l) M(dy) 
im 

^^e^ffCw) (1 _ g-^w) M{ds)^ Ti{x) 

I 

where in the last step we have considered (j2.4p with drift term 6 = 0. Indeed, we are 
considering subordinators with no drift. Also, we have that 

Y^dte-'^^^^'^Tiix) = ^e-**(^') {-^{fii))Ti{x) < oo. (2.23) 
I I 

From this, we arrive at 

{dt - B%)rjf{x) = <^Udt - ]D)^,)e-**(^')y;^(x) = (2.24) 

Im 

end equation ()1.28p is satisfied. 

For the completeness of the set {Y/m : / > 0, m = —/,...,+/}, the series (|1.29p 
converges in the sense that, Vi, 



lim E 

L— >oo 



Lt{^) -Y,e-'^^^'^^Ti{x)\ \{dx] 



Indeed, if we write 

N N 



1=0 1=0 \m\<l 



for a fixed t > 0, then we obtain that 



L 

2i*(w) 



(=0 \m\<l 

=E\\vff-j2 5^(2/ + l)C;e-2**('^') 

1=0 \m\<l 

where 

\\vff=[ \vnx)\'X{dx) and E||r?f f = A(S?)E[ryf (x)]^. 
Js'i 
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Prom p.23p . by recalling that A(Sf) = An and taking the limit for L ^ oo, we get that 

lim E\\i]f - Sif = 0. 

L—^oo 

Furthermore, from ()2.2ip . 

/ Im 

= ^ aimPt {t^l)Ylm{x - Xq) 
Im 

= Y,P^''if^i)Ti{x - xo) = ^e-**(^')rKx - xo). 
I I 

The proof is completed. □ 

Proof of Corollary m We use formula (j2.5p . We assume that ()1.36p holds true and, from 
([222]), we get that 



-(-As2)%(x) =,^^^^j^ {P.,m{x)-m{x) 



ds 



=Y,p^''i^^l)i- fornix) 
I 

=Y,dtP^{pm{x). 



I 



We also observe that P*(/u) = 1 and therefore we get riQ{x) as in the claim. □ 

Now we prove our last result. 
Proof. Proof of Corollary [2] It follows from the proof of Corollary [TJ Prom (|1.36p we can 



write 

oo oo 



Prom (jl.lSp . by setting ti = t2 = t/2, we get that 

oo 

jtix,y) = Y,^~"''nTi{x)Ti{y)] (2.26) 
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and therefore, 



=E 



1=0 

oo +1 

Yl aimaim'{dt + {-As2T)e-^^'Yim{x)Yim'i- 



1=0 m,m'=—l 

The same reasoning apphes for 



□ 
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